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' Abstract 

<n : 

The Lax pair representation in Fourier space is used to obtain a list of integrable scalar evolutionary equations 
with quadratic nonlinearity. The famous systems of this type such as KdV, intermediate long- wave equation (ILW), 
Camassa-Holm and A. Degasperis systems are represented in this list. Some new systems are obtained as well. The 
generalizations on two-dimensional and discrete systems are discussed. 

An important class of the integrable evolution systems is a class of systems with quadratic nonlinearity. Tradi- 
tionally, the evolution integrable systems are classified by power of dispersion law. We don't assume a concrete form 
' of dispersion law beforehand. Moreover, we consider dispersionless systems also. 
We shall consider the systems of following type 

4ru q {t) =u(q)u q (t) + ^2 w (Pi^P2)u P i{t)u P2 (t) (1) 
Pi+P2=q 

> 

The summation in ([TJ) has a symbolic meaning. One can consider a finite set of values of " momentum" p - this case 
corresponds to the finite-dimensional systems. The other case p 6 corresponds to the one-dimensional evolution 
■ equations. 

In fact, one can apply the inverse Fourier transform (we will omit an imaginary unit) to equation (Q) to obtain the 
CN ' systems in a coordinate space. 

For example, one can substitute cu(k) = fc 3 , w(p, q) — p + q to the equation (Q) and apply the inverse Fourier trans- 
form to obtain the KdV equation. Another example: the intermediate long-wave water (ILW) equation corresponds 
to a choice uj(k) = k 2 ]^ hk , w(p, q) = p + q in equation (|l|). 

One of the most effective tools used in the theory of integrable systems is Lax pair representation. One can try to 
^ . obtain Lax representation of (0) in the following form 

L t = [A, L] , L = a + J2u p u p , A = /? + ^ T/ p u p , [at,0\=O (2) 



X 



where a, (3, U p , V p - some constant operators. It is easy to see that (0) is Lax pair for p) if and only if 

[V p , U q ] + [V q , U p ] = 2w(p, q)U p+q , \p, U p ] + [V p , a] = tu(p)U p (3) 
One can use a following matrix representation for the operators in Fourier space 

a k k> =a(k)5k,k<, Pkk> =/3(k)5k,k', U p k , — l(k,p)6k,k'+ P , V p k , = a(k,p)S k ,k>+ P (4) 
One can substitute (^) to (||) to obtain 

a(k,q)l(k - q,p) + a(k,p)l(k - p,q) - a(k - q,p)l(k,q) - a(k - p,q)l(k,p) = 2w(p,q)l(k,p + q) (5) 

and 

[ Q ( fc ) _ a (k -p)]a{k,p) = \J3(k) - P{k - p) - u(p)]l(k,p) (6) 
It is easy to obtain a particular solution of (||). 

a(k) = l(k, 0), [3{k) = a(k, 0), u(k) = 2w(k, 0). (7) 
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This solution is compatible with (||), but is trivial because it corresponds to a shift u p (t) — > u p (t) + 6 p> o in ([!]). 
We shall consider l(k,p) = case now, then equation (|^) transforms to 

a ( k,q ^ p(k-q) +a ( k ' p ^ pfk^p) ~ a<yk ~ q ' ^ ~ a<yk ~ P ' ^ = 2wi(P ' ^ ^ 

It is natural to differentiate (^|) on p and on q to obtain 

a(k,p) = kA(p) + b(p)+c(k), 2w(p,q) = P A{q) + qA{p) + B{p) + B(q), b{-p) = b(p) + pA{p) + B{p) (9) 

where A(p) = A(—p), B{p) = —B{—p), uo{p) — —bj{—p). 

One can substitute (||) to (ra) to obtain a functional equation 

U(x,y) = U(y,x), U(x,y) = p(x)[xA(x - y) + b(x - y) + c{x)\ (10) 

that dehnes a whole class of integrable dispersionless (ui p = 0) systems with quadratic nonlinearity because one can 
choose a = 0, (3 = in this case. 

One can rewrite (|^) to obtain another functional equation 

(3{x) - (3{y) - uj(x -y) = (a(x) - a{y))U{x, y) (11) 



where U (x, y) is some solution of ( |10[ ) . 

The general solution of ([ll]) defines a class of integrable systems (jl|) with nonzero dispersion (u> p ^ 0). Evidently, 
this solution doesn't exist for all U(x,y) from ([LOh. 



1 Classification of dispersionless systems 

It is possible to obtain the general solution of (|ic|). 

The first step is to apply operators D = d x + d v , D 2 D 3 to equation (|l0|) to obtain three more equations - it is 
important that D f{x — y) = 0. Therefore, equation ( J10| ) transforms to 

W D (xp(x) - yp{y), p(x), p(y), p{x)c{x) - p(y)c(y)) = (12) 

where D— Wronskian is defined by formula 

W D {f x {x,y), f 2 (x,y),..,f N (x,y))=det(Wi), Wf = D^ 1 fi(x,y), i,j = l,2,..,N (13) 

The second step is to apply y — ► x limit to equation (|l2|) to obtain an necessary condition 



R(x) 

c(x) — pi + p 2 x + P3, \ip"(x) + X 2 xp'(x) + X 3 p'(x) + Mp(x) = 0, R'(x)=p(x) (14) 
p(x) 



There are two cases - case A2 = is "trivial" - any solution of (jlj) corresponds to a solution of ([L2| 
p{x) — Ci e hlX + c 2 e h2X or p{x) = (x — a)e hx or p(x) — e hx or p(x) = x or p(x) = 1 



In A 2 ^ case one can substitute xp{x) — — j-(Xip'(x) + (A 3 — X 2 )p(x) + A4_R(x)) to (12), to take the y — > x limit 
as well to obtain the second necessary condition. 

v A p IV (x) + v 2 p"{x) + vxp'(x) + ^ = (15) 

In this case we have only three solutions different from A2 = case: p(x) — (x — a)(x — a — b)(x — a + b) or 
p(x) — (x — a) (x — b) - Degasperis and Camassa-Holm systems respectively and p(x) = 1/x case. 

So we have obtain following classification of dispersionless systems with quadratic nonlinearity: 

1. p(k) = 1 

In this class one can found the hydrodynamic type systems 

u t = mu x , m p =A(p)u p , L kk * =u k -k', A kk * = ^(k + k')A(k-k')u k -k>- (16) 

2. p(k) f 1, A(k) = 

• The system from this class have a general form in a p(x) = e hx case 

u t =mu, m p =B(p)u p , L kk > = e~ hk u k - k ,, A kk , = b(k - k')u k -k', B(p) = b(p){e' lp - 1). (17) 

It is interesting that a choice b(x) — 1+ ^ x corresponds to Hilbert-Hopf equation vnt — L(d x )m 2 , L(k) = 
th^ where (it* = m k (e hh + 1)). 



• The case p{x) = (x — a) e hx corresponds to 

1 e Kk'-k) _ i 

u t = mu, m p = B(p)u p , L kk > = (k - a) e x u k - k ,, A kk > = ( - - - H - — — )u k -k> ■ (18) 

where B(p) = ^(2 - e hv - e~ h P) 

• Two-exponent case p(x) = e hlX + e h2X corresponds to 

5 t n - e' ll2: ) - 6 2 (1 - e' 12 ^) 
u t = mu, m p = B(p)u p , B{jp) = b{-p) -b{jp), b(x) = ^ _ ^ (19) 



where L fcfc - = (aie' ,lfc + a 2 e h2k )- 1 u k - k ', A kk , = (b(k - k 1 ) + c{k))u k - W , c{k) 



a 1 e h i k +a 2 e h 2 k 



Note that the systems ( |18[ ) and (19) have a second Lax pair representation (|l 
3. P (fc) ^ 1, A(k) ? 

• First example in this class is a system with p(x) = x 2 — 1/4 : Lax pair has a form 

= k 2 -l/A Uk - kl ' Akk ' = \ { \ k {k ik k'Y + fc2-l/4 )Ufe - 
that corresponds to Camassa-Holm H equation 



u t = 2f x u + fu x , u=-(f xx -f) (20) 



• Second example is a system with p(x) = x{x 2 — 1) : Lax pair in this case has a form 



1 . k-2k' fc 

Lfcfc '-jfe(fc2-i) ^ fcfc ' - 4 _ (fe - to 2 + fe^T^*-*'- 

that corresponds to A. Degasperis system 

u t = 3f x u + fu x , u=-(f xx — f) (21) 

The other cases are 

Systems with /?(a;) = x and A kk > = k' A(k — k') have a form 

u t = d x (mu), m p = A(p)u p (22) 
Systems with p{x) = - and A kk ' = kA(k — k') have a form 

Ut — mu x — m x u, m p — A(p)u p (23) 



2 The systems with dispersion 

One can try to solve equation (|TT|) starting from the general solution of (JToh to obtain a list of systems with dispersion. 



One can verify that there are no solutions of (|llj) for systems (|20 
obtain a generalization of Camassa-Holm (^o|) and A. Degasperis 
is valid for the systems ( |l8|) and (|l9| ) as well. 

The rest of possibilities are systems with p(x) — 1, p(x) = e hx or p(x 



and (|21|). So the statement is: It is impossible to 
2l| ) on the systems with dispersion. This statement 



1. In a case p{x) = 1 it is easy to obtain U(x,y) = ■k(x + y)A(x — y) +d(x — y), where A(x), d(x) are even functions. 
One can substitute given U(x, y) to (pr[), apply the operator D to equation to obtain a necessary condition 



W D (a(x) - a(y), (x + y)(a(x) - a(y)),/3(x) - f3(y), 1) = (24) 



The method of solution of equation ( |24|) is similar to that of solution of fll2|): Limit x — * y gives 13' {x) = cxa'{x). 
Substitution (3(x) to (|24|) and taking x — > y limit give a necessary condition a" (x) + cixa'(x) + c^Oi'ix) + C3 = 0. 
The only possible choice to fulfill (|i|) is a(x) = x 2 , (3{x) = 4x 3 (one can choose c = 6). 

It is easy to obtain u>(x) — x 3 in this case, that corresponds to KdV equation u t — u xxx + 6uu x with 



k 2 5 k k' + u k - k >, A kk t = 4fc 3 <5 fc fc' + 3(fc + k')u k - k r. (25) 



2. In p(x) = e hx case one have U(x,y) — e^ x+v ^d{x — y), where d(x) is even function. 
The necessary condition in this case 

W D ((a(x) - a(y))e^+y\(3(x) - f3(y), 1) = (26) 

One have (3'(x) — Cia'(x)e hx +c 2 . The only non-trivial solution of ( ]26| ) a(x) = x e~ hx , (3(x) = x 2 gives interme- 
diate long-wave (ILW) equation 

u t = T{u xx ) - 2uu x , r(p)=coth(^) (27) 



with Lax pair 

k — k' 

L k k' = ke- hk 5 kk , + e~ hk u k - k ,, A kk , = k 2 5 kk , + 2 — h{k _ kl) u k ^ k , (28) 
3. In p(x) = x case we just obtain alternative Lax pair for KdV equation: 



k p 1 k jy^ ^ 

L = 5 kk ' {- + -) + -u k - k ,, A = S kk '(—-fj,k) + k'u k -k', L t = [A,L] u p . t = —u p +p22u q u p - q . (29) 



3 Generalizations 

It is possible to extend our approach to high-dimensional and multi-component systems. Consider two-dimensional 
scalar systems with quadratic nonlinearity. Unlike an one-dimensional case, the "wave" variables k,p and q in the 
functional equations (|^) and (^) are vectors with two components. To obtain the functional equations in this case one 
can use following trick: One can write k = (k x , k y ) = (k, A), where k is an argument of the functional equations and 
A is parameter. For example, one can rewrite equation ( [l0|) as follows 

p(x, A) [xA(x - y, A - p) + XA 2 (x - y, A - p) + b(x - y, A - p) + c(x, A)] = . . 

P(y, /") [yA(y -x,p-X) + pA 2 (y-x,p-\) + b(y- x, p - A) + c(y, p)] 



The method of solution of ( p0[ ) is analogous to one-dimensional ones. One can use notations p(x) = p(x, A), c(x) — 
c(x, A), p(y) = p(y, p), c(y) = c(y, p) to obtain 

W D {xp(x) - yp(y), p{x), p(y), p(x)c(x) - p{y)c{y)) = (31) 



in this case (compare with (|12|)). 

Next step is to take p — ► A limit in equation (^lj) to obtain equation ( |l2| ) exactly! Note that one can't apply this 
limit to functional equation ( |30| ) because possible divergence in the functions A, A 2 , b. 

Then the only possibility to obtain two-dimensional systems with quadratic nonlinearity is to start from one- 
dimensional systems and to try to extend some solutions of ( |l0| ) to (p0|). We have no full classification in two- 
dimensional case at the moment, but a couple of examples are in agree with our approach: 

First example is a famous Kadomtsev-Petvishvili (KP) equation 

u t = u xxx + 36» _1 (M y ) + 6uu x , O^ 1 = dyd' 1 

has a Lax pair in Fourier space 

k — k' 

hk' = (k y + k 2 x )5 kk > + u k -k', akk' = 4fcf4fe' + 3(k x + k' x - -^-)u k - k > (32) 

K x K x 

that is an agreement with the general form of Lax pair (^) in our method. 
Second example is Veselov-Novikov (VN) equation 

u t = u xxx + u vvv + u9(u x ) + U x 6(u) + uQ~ l (uy) + u y 9^ 1 (u), 9 = d x dy X (33) 



One can consider VN equation (|33j) as a sum of two symmetries u t = u t + + u t - where u t + = u xxx + u9(u x ) + u x 9(u) 
and u t - = Uyyy + u6~ l (uy) + u y 9~ 1 (u). 

Lax pair in Fourier space for " +" symmetry has a form 

1 k — k' 

L kk ' = 3fc y 4fc' + — Ufe-fc', A kk , = kl5 kk > + k' x -^ — jjUk-k' (34) 

K x Ky Ky 



One can introduce new operators: 

Li = L, L 2 = L[x <-> y], A x = A, A 2 = A[x <-> y], L itti — [L l ,A l ], L = L x + L 2 , A = A x + A 2 , (35) 
It is easy to obtain the Manakov's triada representation for Veselov-Novikov equation (|3£ 

L t = [A, L] + BL, B = A{ 1] + A [2) -A Y -A 2 (36) 

A case of multi-component systems corresponds to a choice fc — > (fc,n) where k 6 5ft and n = 1,2, ..,JV. A 
classification of those systems is a subject of forthcoming papers, but in a case of two-component systems we have two 
interesting classes: 

First example is so called Toda-type systems 

z q ,t = ~u{q)z q + y(q)p g + ^r (1) (s,<7 - s )z s z q ^ s , p q>t = u>(q)p q + ^T (2) (s,q - s)z s p q ^ s 

s s 

with Lax pair 

L = a + J2u q z q + Y,W q Pq A = p + J2V 9 z g , [a,0\=O (37) 

9 9 9 

Note that Toda system has Lax pair representation in coordinate space 

L t = [A,L], L = d + z-X+pd~ 1 , A = d 2 + 2d z <^ z t = d\-z + z 2 + 2p], p t = d\p' + 2pz}. (38) 

This dynamical system can be transformed to nonlinear Schrodinger equation (p = i/)*ip, z = —d x log(ip)) and 
Backlund transformation of this system have a form of Toda lattice (see Q for more details). It is easy to see that 
Fourier transform of (|38| ) has a form (p7|) . 

The other example is a system with 3- wave interaction 



i a p>t = uj(p)a p + W {1 \p, q)a p+q a* + £ (p, q)a q a 



P— 9 



>(p)a; + £ VF (1) (p, <z)a; +g a, + £ W>( 2 ) (p, <z)a* a ;_ g (39) 



The problem is to find all possible functions and u> in case when system ( |39"| ) is integrable. This question 

was discussed in || by the use of perturbative approach. 
Lax pair in this case has a form 

L = a + Y / U p a p + Y / (U p ) + a;, A = + £ V*a p + ^(^ p ) + a;, [«, >S] = 0. (40) 

p>0 p>0 p>0 p>0 

One can derive a system of functional equations as in scalar systems case, but a general solution of these equations is 
unknown as yet. It is interesting that system ( J3S| ) is equivalent to scalar system ([!]) in a "real" case: 

u p = 0(p)a p + 6(-p)a*_ p , u* = W {1 \p, q) = w{q,p - q), W {2 \p,q) = w{p + q,-q). 

In summary, we have use a Fourier representation of Lax pairs to classify scalar integrable systems with quadratic 
nonlinearity. The similar (but not the same) approach is so called symbolic method. This method was applied to 
theory of integrable systems by Gelfand and Dickii || and was improved in the recent papers W],|§]. The common 
feature and main advantage of both methods is that all coefficients of equations in Fourier space are the functions 
(not operators!) of wave numbers k. Because this fact, all necessary conditions can be formulated as the functional 
equations. The irony of it is that in some cases it is better to solve functional equations rather then differential ones. 

This work was supported by INTAS No. 99-01782 and Russian foundation for Basic Research No. 01-01-00874. 
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